Astronomy & Astrophysics manuscript no. Franz ' Kirsten 


©ESO 2013 


March 2,2013 





No evidence for a central IMBH in M15 

F. Kirsten 12 '* and W.H.T. Vlemmings 3 

1 Argelander Institut fur Astronomie (AlfA), University of Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 

2 Max Planck Institut fiir Radioastronomie (MPIfR), Auf dem Hiigel 69, D-53121 Bonn, Germany 

3 Department of Earth and Space Sciences, Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, 
Sweden 



ABSTRACT 

Intermediate mass black holes (IMBHs) with expected masses M. ~ 10 4 M are thought to bridge the gap between stellar mass black 
holes (M. a 3 - 100 M e ) and supermassive black holes found at the centre of galaxies (M. > 10 6 M e ). Until today, no IMBH has 
been confirmed observationally. The most promising objects to host an IMBH as their central mass are globular clusters. Here, we 
present high sensitivity multi-epoch 1 .6 GHz very long baseline interferometry observations of the globular cluster M 15 that has been 
suggested to host an IMBH. Assuming the IMBH to be accreting matter from its surrounding we expect to detect it as a point source 
moving with the global motion of the cluster. However, we do not detect any such object within a radius of 6000 AU of the cluster 
centre in any of the five observations spread over more than one year. This rules out any variability of the putative IMBH on the time 
scale of one to two months. To get the most stringent upper limit for the flux density of the putative IMBH we concatenate the data of 
all five epochs. In this data we measure a 3c upper flux limit of 10 jijy for a central source. We employ the fundamental plane of black 
hole activity to estimate the mass of the central IMBH candidate. Based on previous X-ray observations of M15 our measurements 
indicate a 3tr upper mass limit of « 500 M . 
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1. Introduction 

Super-massive black holes (SMBHs) with masses M. a 10 6 ~ 9 
solar masses (M Q ) are known to exist at the centre of galax- 
ies (Kormendy & Richstone 1995). At the same time, obser- 
vations of high-mass X-ray binaries led to the conclusion that 
stellar-mass black holes with M. a 100 M Q must form as well 
(McClintock & Remillard 2006; Ozel et al. 2010). The exis- 
tence of black holes (BHs) bridging the gap between these two 
extremes in mass, the so-called intermediate mass black holes 
(IMBHs, M. * 10 4 M@), is still under debate. 

One possible candidate for IMBHs are ultraluminous X-ray 
sources (ULX, Colbert & Mushotzky 1999) appearing to accrete 
matter at super-eddington rates. The interpretation of ULX- 
properties as being characteristic for IMBHs, however, is still in 
discussion (e.g. Berghea et al. 2008, Zampieri & Roberts 2009). 
Nevertheless, the recent discovery of the hyper-luminous X-ray 
source ESO 243-49 HLX-1 by Farrell et al. (2009) adds evidence 
that ULXs might host IMBHs. 

Based on the black hole mass-stellar velocity dispersion re- 
lation (M. - <x) established for galaxies (Ferrarese & Merritt 
2000; Gebhardt et al. 2000), an object like an IMBH can be 
expected to reside at the core of globular clusters (GCs). Even 
though recent work by Vesperini et al. (2010) shows that an 
IMBH can be formed in a GC environment, no conclusive ob- 
servational evidence for the existence of IMBHs has been found 
to date. 

Apart from, e.g, the globular clusters to Cen (van der Marel 
& Anderson 2010), Gl in M31 (Gebhardt et al. 2002; Ulvestad 
et al. 2007), and 47 Tuc (Lu & Kong 201 1), M15 has been one 
of the most promising GC-candidates to host an IMBH for a 
long time. Stellar surface density profiles of M15 reveal a steady 



increase towards the center of the cluster indicating a state of 
advanced core-collapse (Djorgovski & King 1986). To explain 
the central brightness peak Newell et al. (1976) suggested an 
800 M black hole to reside at the core of the cluster. The in- 
crease of the stellar velocity dispersion towards the centre of 
M15 as measured by, e.g, Gerssen et al. (2003), supports the no- 
tion of a high central mass concentration in the cluster. In fact, 
dynamical models based on line-of-sight velocities and proper 
motions infer a mass of 3400 M within the central 1 arcsec- 
ond (= 0.05 pc at the distance of 10.3 + 0.4 kpc, van den Bosch 
et al. 2006). The nature of this mass concentration is unknown. 
Gerssen et al. (2003) invoke the existence of an IMBH with 



a mass of M. 



1700 H 



1700 



M G to explain their observations. 
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Similarly to Illingworth & King (1977), Baumgardt et al. (2003) 
and McNamara et al. (2003) challenge this interpretation, based 
on N-body simulations excluding an IMBH. Instead, these au- 
thors claim that the observational data can be explained equally 
well by assuming a collection of neutron stars to exist at the 
core of M15. The latest simulations of that kind require a total 
of 1600 neutron stars to fit the velocity-dispersion profile men- 
tioned above (Murphy et al. 201 1). 

Complementary to the indirect (non-)evidence using kine- 
matic studies based on optical observations, X-ray and radio ob- 
servations aim at directly detecting such an object. The funda- 
mental plane of black hole activity (FP) as determined for active 
galactic nuclei (AGN) (Merloni et al. 2003; Falcke et al. 2004), 
relates black-hole mass, X-ray and radio luminosity. Now, as- 
suming the same physical processes powering AGN-emission 
to also be characteristic for IMBHs, Maccarone (2004) predict 
a black hole mass of 400 M Q for the IMBH in M15. Their re- 
sult is based on an estimate of the cluster's total mass, Mac, 
from its absolute V-magnitude, My, and on a model by Miller & 
Hamilton 2002 stating that M. * 10" 3 M GC - Bash et al. (2008), 
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on the other hand, perform a survey of Ml 5 at 8.6 GHz using 
the Very Large Array that reaches a noise level of 8.5 yuJy/Beam. 
Using the FP they predict a flux density of 10 3 - 10 5 pJy for a 
putative IMBH with a spectral index a = -0.7. They detect no 
central source at a 3cr upper flux limit of 25 juJy. Similary, Cseh 
et al. (2010) tried to detect an IMBH at the center of the globular 
cluster NGC 6388. The 3<x noise level (81 piy) of their observa- 
tions with the Australia Telescope Compact Array allowed the 
authors to constrain the mass of the possible IMBH to be lower 
than * 1500 M Q . 

In this paper we discuss multi-epoch observations of M15 al- 
most three times as sensitive as those of Bash et al. (2008). The 
high angular resolution of our data allows us to disentangle any 
possible background sources from objects belonging to the clus- 
ter. Furthermore, the long time line spanning 15 months allows 
us to, in principle, detect the proper motion of a possible central 
radio source moving with the global motion of the cluster ex- 
pected to be on the order of -1.0 + 0.4 and -3.6 + 0.8 mas/yr in 
right ascension (RA) and declination (Dec), respectively (Jacoby 
et al. 2006). 

2. Observations 

We observed M15 five times in a global VLBI campaign that 
was spread over a time period of more than one year. The 
observations included in this analysis were conducted on 11 
November 2009, 7 March 2010, 5 June 2010, 2 November 
2010, and 27 February 2011. The array we employed con- 
sisted of eight European VLBI Network (EVN) antennas (Jodrell 
Bank, Onsala, Westerbork, Effelsberg, Noto, Medicina, Toruh, 
Arecibo) and the Greenbank Telescope (GBT). We observed at a 
central frequency of 1.6 GHz and the data was recorded at 1024 
Mbps. Accounting for the different receiver systems at the indi- 
vidual telescopes our total bandwidth amounts to 230 MHz on 
average. The correlation was done at the EVN-MklV correlator 
(Schilizzi et al. 2001) at the Joint Institute for VLBI in Europe 
(JIVE). 

The longest baselines in east-west (north-south) direction ex- 
tending over 7500 (2000) km allow for a resolution of 2.2 X 6.3 
mas. The largest dishes of the array (Arecibo, Effelsberg, and 
the GBT) ensure a maximum sensitivity of approximately 4 
/Jy/Beam. 

The observing schedule lasted six hours in total, 3.6 hr 
of which were spent on the target cluster Ml 5. The quasar 
J2139+1423 (located « 3.17° to the north-east of the pointing 
center) served as phase calibrator and the blazar 3C454.3 was 
used for bandpass calibration. Arecibo-data is available for 75 
(50) min in epochs 1 and 5 (epoch 3). Unfortunately, epochs 2 
and 4 lack any Arecibo data which is why the sensitivity and 
astrometric precision of these two datasets is lower by about a 
factor of two. 

Aiming to detect compact radio sources close to the core 
of M15 we map out the entire central region within 2' in only 
one pointing. For this project, however, only the very central re- 
gion (the central 16") correlated at RA = 21 h 29 s 58*3120, Dec = 
12°10'02'.'679 (J2000 equinox) is of interest. The entire dataset 
will be described in a forthcoming paper. 

3. Data reduction 

After correlation, all data is reduced, calibrated and imaged us- 
ing the NRAO Astronomical Image Processing System (AIPS 1 ). 

1 http://www.aips.nrao.edu/ 



A priory calibration tables including system temperature and 
gain curve corrections as well as a flag table containing in- 
formation about band edges and off-source times are provided 
by the EVN pipeline 2 . We apply these to the dataset as given. 
Parallactic angle corrections are determined with the AIPS task 
CLCOR and first ionospheric corrections are computed running 
TECOR with the total electron content (TEC) maps published 
by the Center for Orbit Determination in Europe^. Even though 
these maps are quite crude in angular resolution (about 5° x2. 5°) 
they have shown to be of use reducing the scatter in phase delay 
by a factor of 2-5 (Walker & Chatterjee 1999). Next, we identify 
and flag radio frequency interference (RFI) for all antennas and 
sub bands. 

The bandpass calibration is done running BPASS on the data 
for 3C454.3 and yields phase and amplitude gain factors for all 
8x128 (512, epoch 1) channels for all antennas. We align phases 
in between IFs by performing a manual fringe correction running 
FR1NG on 3C454.3 on a sub-interval of about 30 s of observa- 
tion. 

At this point, we combine all correction tables obtained so 
far and apply it to the data of the phase calibrator. We fringe fit 
this dataset including data over the entire time range. We solve 
for phase delays and phase rates simultaneously using solution 
intervals of 1.5 min. 

The fringe solutions in conjunction with all calibration so- 
lutions found earlier are then applied to the M15 data. In or- 
der to eliminate any residual phase delays and amplitude errors 
caused by the atmosphere and the ionosphere we take advantage 
of the strong unclassified source SI (Johnston et al. 1991) lo- 
cated about 94" to the west of the cluster centre and use it for 
in-beam calibration. To speed up the self-calibration process we 
average the visibilities both in the time- and frequency domain 
to 2 s integration time and 64 channels per IF. 

Finally, we image the self-calibrated data running IMAGR 
employing natural weighting to ensure maximal sensitivity. In 
order to account for possible inaccuracies in the assumed cluster 
centre we produce an image that has an angular size of roughly 
16" x 16" (= 0.8 x 0.8 pc). 

4. Analysis and results 

In all five epochs we first produce a noise map of the image by 
smoothing it with a kernel that has a size of 1024 x 1024 pix- 
els. Based on this rms map, the AIPS source detection algorithm 
SAD then searches for objects down to a signal-to-noise ratio 
of 3. The coordinates of the possible sources detected in this 
fashion in all five observation epochs are then cross-correlated. 
Cross-correlation is performed allowing for a maximal posi- 
tional shift of 15 mas in between epochs. This corresponds to 
roughly four times the maximal beam width in right ascension. 
No match can be found relating all five data sets. 

Finally, we also inspect the images manually. Figure 1 dis- 
plays contour plots of the dirty images of the central region 
of M15 for all five epochs. All tiles are centred on the clus- 
ter core as published by Goldsbury et al. (2010) at coordinates 
RA = 21 h 29 s 58 s .330, Dec = 12°10'01'/200 (accurate to within 
0.2"). We do not detect a significant signal in any of the epochs. 
Accounting for the different sensitivity limits varying between 
4.3 juJy/Beam in epoch 3 and 11.5 juJy/Beam in epoch 4, we 
can put upper 3cr limits between 13 and 35 //Jy/Beam on the 
flux density of a possible central object. Figure 1 also displays 

2 http : //www . evlbi . org/pipeline/user_expts . html 

3 ftp://ftp.unibe.ch/aiub/CODE/ 



2 



F. Kirsten and W.H.T. Vlemmings: No evidence for a central IMBH in M15 



1.0 FF 



0.6 
0.4 



^ 1.0 FF 



Epoch 1 




- 0.8 - 

- 0.6 - ' 

- 0.4 

- 0.2 - 
,- 0.0 

-0.2 



-i.o kt 



■'Epoch 1 




"J " • 



1.0 FF 



g o.o 
< 



Epoch 3 




E*podj4 



-1.0 fc_ 
1.0 




0.6 4'' 
0.4 
0.2 - 

o 

LU 

V) 

o oo 
< 

-0.2 



Epoch 5 ; 



n=i i.o 



FF~ 




0.0 

ARC SEC 



-0.2 
-0.4 
-0.6 
-0.8 



Coiriposite 




0.0 

ARC SEC 



Fig. 1: Contour plots of all five individual observations and the composite image combining the data of all epochs. The tiles are 
centred on the assumed core of M15 at coordinates RA = 21 h 29 s 58 ! 5330, Dec = 12 o 10'01'.'200. The circle indicates the 3cr error of 
the position of the core as determined by Goldsbury et al. (2010). Adopting a distance of 10.3 kpc to M15, each tile has a physical 
dimension of about (20000 AU) 2 . Contours are (-5, -3, 3, 5) times the rms in each individual epoch (4.7, 8.9, 4.3, 1 1.5, 5.8 yuJy). 
The noise level of the composite image is 3.3 fiJy. 
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a dirty image of the data concatenated over all five epochs. The 
rms of this deconvolved image is 3.3 juJy/Beam which translates 
to a 3cr upper flux limit of 10 //Jy/Beam. 



5. Discussion 

5. 1 . Radio flux limits and variability 

The five individual epochs have rather different sensitivity limits 
that correspond to a 3cr upper flux density limit of 14.1, 26.7, 
12.9, 34.5, 17.4 /iJy/beam for epochs 1 to 5. We did not detect 
a central source in any of the observations that were conducted 
at regular intervals of about three months spanning a time range 
of 15 months. Provided the IMBH-candiate is of transient nature 
that is in its 'on'-state for longer than one or two months, the 
probability to have missed it in all five observations is negligible. 
Therefore, we rule out any variability of the central object on 
these time scales. Accordingly, we can assume a steady state 
emission model in which case the noise level of the concatenated 
data yields the most stringent 3cr upper flux limit of w 10 fiJy. 
This flux limit is a factor of 2.5 lower than that from Bash et al. 
(2008) and we will adopt it throughout the following analysis. 
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5.2. Mass limits from X-ray observations 

Similarly to the first versions of the FP from Merloni et al. (2003) 
and Falcke et al. (2004), the FP derived by Kording et al. (2006) 
uses a sample including both X-ray binaries (XRBs) and AGN. 
Thus, both versions span several orders of magnitude in BH 
mass and should also be applicable for the intermediate mass 
range of IMBHs. The relation found by Kording et al. (2006) 
for XRBs and low-luminosity radiatively inefficient AGN has 
the lowest intrinsic scatter cr int = 0.12 dex (^30 percent) and we 
will use it in the following. In terms of black hole mass M„ radio 
luminosity L«, and X-ray luminosity Lx their FP-relation reads 

logM. = 1.55 \ogL R - 0.98 logL x - 9.95 . 

Ml 5 is known to host two strong X-ray sources, AC211 
(Giacconi et al. 1974; Clark et al. 1975) and M15 X2 (White & 
Angelini 2001), close to the core of the cluster. Both have been 
classified as low mass XRBs. The strong X-ray emission of these 
two objects makes it difficult to detect the expectedly faint emis- 
sion of a central IMBH. Nevertheless, Ho et al. (2003) put an 
upper limit of Lx — 5.6 x 10 32 erg/s on the X-ray luminosity of 
the putative IMBH. Hannikainen et al. (2005) manage to detect 
a faint source (Lx = 3.3 x 10 32 erg/s) close to core of M15 which 
they attribute to a dwarf nova. In their paper, Hannikainen et al. 
(2005) report a flux detection limit of 2 x 10~ 15 erg/s/cm (0.5- 
2 keV) which translates to an X-ray luminosity of 2.54 x 10 31 
erg/s at the distance of M15. Converting our measured radio flux 
to a 5 GHz radio luminosity Lr < 6.34 x 10 27 erg/s (assuming a 
flat radio spectrum) and inserting it together with the upper limit 
for the X-ray luminosity from Hannikainen et al. (2005) the FP 
yields a mass limit of 232 + 93 M G . This translates to a very con- 
servative 3cr upper mass limit of 5 1 1 M Q for the putative IMBH. 
Our observations thus decreased the upper mass limit by a fac- 
tor of four compared to that indicated by earlier observations by 
Bash et al. (2008). 

5.3. Mass limits from accretion models 

The lack of a central X-ray source in M15 is expected because 
the accretion process is most likely radiatively inefficient (Ho 



Fig. 2: Black hole mass as a function of X-ray luminosity for 
our measured radio luminosity (Lr = 6.34 x 10 27 erg/s) as pre- 
dicted by the FP from Kording et al. (2006) (solid lines) and 
Plotkin et al. (2012) (dashed lines). The gray shaded areas in- 
dicate the lcr uncertainties of the relations. The horizontal lines 
mark the predicted masses for different assumptions on the ra- 
diative efficiency 77 and the accretion efficiency e. The vertical 
dotted line marks the upper limit for the X-ray luminosity from 
Hannikainen et al. (2005). 

et al. 2003). To constrain the mass of the IMBH-candidate we 
estimate the X-ray luminosity based on the (unknown) accretion 
rate M: 

L x = r]e c 2 M , 

where r\ and e are the radiation efficiency and the accretion ef- 
ficiency, respectively. The accretion process of radiatively in- 
efficient black hole binaries (e.g. Esin et al. 1997) as well as 
those of quiescent SMBHs (e.g. Yuan et al. 2003) is well de- 
scribed by advection-dominated accretion flow models (ADAF, 
Narayan & Yi 1994). Such BHs undergo quasi-spherical accre- 
tion and, thus, following Maccarone (2004) we further assume 
that mass accretion of the putative IMBH in M15 can be de- 
scribed by the Bondi-Hoyle-Lyttleton (BHL) formalism (Hoyle 
& Lyttleton 1941; Bondi & Hoyle 1944; Ho et al. 2003) 

^ = 32X 1017 (200^) 2 (cl2^) (wkJ 15 g S l ' 

Here, « and T are the gas density and temperature in the GC, 
respectively. We adopt n = 0.2 cirT 3 from Freire et al. (2001) 
and use a typical GC gas temperature of T = 10 4 K. Inserting 
the above equations into the FP, and solving for M, yields 

2.96 logM. = 1.55 logL s - 0.98 log(77 * e) - 41.17 . 

Mass estimates from this relation depend crucially on the choice 
of J] and e. Observations of the IMBH-candiate in Gl in M31 
suggest that the radiative efficiency 77 < 0.01 for this source 
(Ulvestad et al. 2007). This is consistent with the aforemen- 
tioned ADAF models that are only valid for 77 <k 0.1 (Narayan 
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& McClintock 2008). In order to cover as large a parameter 
space as possible we choose 77 = 0.1 as conservative upper 
limit for the radiative efficiency. As lower limit we take the es- 
timate of Ho et al. (2003) who state that 77 can be as low as 
10~ 4 for the IMBH candidate in Ml 5. It is not clear at all at 
what fraction e of the Bondi-rate the central source is accret- 
ing. Typical values are in the range e = [10 3 , 0.1] (Maccarone 
2004; Ulvestad et al. 2007; Cseh et al. 2010) which we adopt for 
our analysis. Consequently, we explore a parameter space cov- 
ering 77 * e = [0.01, 10~ 7 ]. Based on the upper X-ray luminosity 
as measured by Hannikainen et al. (2005) we can exclude val- 
ues of 77 * e > 10~ 3 because they would yield stronger X-ray 
emission than what is observed. The mass limit for a value of 
77 * e = 10~ 4 is indicated in Fig. 2 and would yield a black hole 
mass M. -94 + 37 M . For 77 * e = 10~ 5 we find M. = 202 ± 80 
M Q which would indicate a source that is accreting at ten per- 
cent of the Bondi-rate with a radiative efficiency 4 77 = 10~ 4 . The 
combination of both, very inefficient radiation (77 = 10~ 4 ) and 
inefficient accretion (e = 10~ 3 ) is also indicated in Fig. 2 with 
the mass estimate of M. = 927 + 371 M Q . 

5.4. Mass estimates from other FP-relations 

Mass estimates computed here with the FP-relation from 
Kording et al. (2006) agree well within the errors with those esti- 
mates computed from the relations found by, e.g., Merloni et al. 
(2003) and Giiltekin et al. (2009). We note, however, a signifi- 
cant offset towards higher masses when using the latest version 
of the FP published by Plotkin et al. (2012). Their relation reads 

logM. = 1.64 logL« - 1.13 \ogL x - 6.89 . 

In Fig. 2 the predicted BH mass from this relation for our mea- 
sured radio flux density is indicated by the dashed lines. With 
this relation the limit for the X-ray luminosity results in a BH 
mass M. = 1654 + 661 M Q . However, this mass limit can only 
be explained by the least efficient accretion limit of 77 * e = 10~ 7 . 
All higher values of 77*6 would yield an X-ray luminosity beyond 
the one observed and can be excluded within this FP-relation. 



6. Conclusions 

Using our multi-epoch high sensitivity observations of M15 we 
were able to put an upper limit of 10 //Jy on the 1.6 GHz radio 
flux density of a central source in this globular cluster. Assuming 
that the central mass concentration is a black hole we employed 
the FP as derived by Kording et al. (2006) to put constraints on 
the mass of the source. The lack of a detection of an object in 
the X-ray observations coinciding with the assumed cluster cen- 
ter allowed us to use an upper limit for the X-ray luminosity of 
the putative IMBH which yields a black hole mass M, = 232+93 
M Q . Such a mass estimate is in agreement with a source accret- 
ing matter from its surrounding at up to ten percent of the Bondi- 
rate with a radiative efficiency as low as 10~ 4 which is consistent 
with results for the IMBH-candidate in Gl in M31 (Ulvestad 
et al. 2007). Even the 3<x upper limit of a 500 M Q , however, 
is still a factor of seven lower than that required to explain the 
dynamics of the cluster (van den Bosch et al. 2006). If we em- 
ploy the FP from Plotkin et al. (2012) the 3cr upper mass limit 
is increased to a 3600 M Q which would agree with the required 

4 Obviously, other combinations of (77, e) would also be suitable but 
we give this one as a limiting case that has been discussed by other 
authors (e.g. Ho et al. 2003, Cseh et al. 2010). 



dynamical mass. However, such a mass can only be explained if 
the central object is accreting at 0. 1 percent of the Bondi rate and 
if only 0.01 percent of radiation is escaping from the accretion 
region. Furthermore, this mass estimate disagrees by a factor of 
seven with results obtained from most other FP-relations and we 
adopt the mass limit of M. = 232 + 93 M e . 

Therefore, we conclude that Ml 5 most likely does not con- 
tain an IMBH but that its central region probably hosts a collec- 
tion of dark remants such as neutron stars as proposed by, e.g., 
Baumgardt et al. (2003); McNamara et al. (2003), and Murphy 
et al. (201 1) instead. The reason for not detecting these pulsars 
is most probably their expected low flux density of w 2 /jjy (Sun 
et al. 2002). 

During the final stages of the refereeing process of this paper, 
Strader et al. (2012) published a mass estimate for the IMBH in 
M15 that is in good agreement with the estimate derived here. 
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